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Abstract 



We exhibit the strong influence on light propagation of the finite size in 
photonic band-gap material. We show that light emission can be controlled 
by the symmetry group of the boundary of the finite device. These results 
lead simply to important practical applications. 

Pacs: 42.10 Qs, 11.80-j, 42.82.-m, 03.50.De, 15.40.Mg 

Photonic crystals are expected to permit the realization of devices for integrated optics or 
laser microcavities The today technology now authorizes the making of bidimensional 

photonic crystals in the optical range 0. From the theoretical point of view, there are 
two ways for characterizing PBG materials. The first way is the one commonly used in 
solid state physics where pseudo-periodic boundary conditions are imposed leading to Bloch 
waves theory. This is a powerful tool for the computation of the band structure of the 
material 0. However, this theory cannot be used when dealing with the scattering of an 
electromagnetic wave by a finite device (which is the actual experimental situation) 



In that case, boundary effects cannot be skipped and must be taken into account. In this 
letter, we address the problem of the effect on light propagation of both the symmetry group 
of the lattice and the symmetry group of the boundary of a finite piece of the material. The 
general study of the symmetry of the lattice modes has already been addressed in order to 



understand the possible degeneracies and the existence of full band gaps fii]HT5||. Here, we 
show that the finite size of the device strongly modifies the behaviour of the electromagnetic 
field, and that the symmetry of the boundary is a crucial parameter for the control of light 
emission from PBG devices. 

We deal with a bidimensional photonic crystal that is made of a finite collection of parallel 
dielectric rods of circular section. The rods are settled according to a lattice with some 
symmetry group Gy (Y denotes the elementary cell) of the plane. The relative permittivity 
of the lattice is a F-periodic function e r (x, y). The rods are contained in a domain Q, having 
a symmetry group Gq. As we deal with objects that are embebbed in an affme euclidean 
space, both groups must be given as subgoups of the group of plane isometries 0(2), in a 
canonical oriented basis. This is due to the fact that the use of the abstract groups does not 
permit to distinguish between two isomorphic realizations. Indeed, using the abstract groups 
and unspecified representations of degree 2, we could not, for instance, distinguish between 
two squares C 1)2 deduced from one another through a rotation r of angle 7r/4: denoting Ti 
a realization of -D4 as an invariance group of C\, then r 2 = rT\r~ l is a representation of C 2 
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and therefore the equivalence class of Ti is not sufficient to compare the two figures and a 
canonical basis has to be precised. 

Denoting 1q the characteristic function of Q (which is equal to 1 inside Q and elsewhere), 
the relative permittivity is given by Eq (x, y) = 1 + ln {e r ( x i y) ~ !)• Assuming a s-polarized 
incident field (electric field parallel to the rods), the total electric field verifies the d'Alembert 
equation: 

c- 2 d tt E z = e^AE z (1) 

Our aim is to study the invariance of this equation under the action of the various groups 
characterizing the geometry of the problem. For arbitrary cross-sections of the fibers, we 
should also introduce their symmetry group (for instance the group -D4 for square rods). The 
choice of circular rods simplifies the study in that their symmetry group O (2) contains both 
Gy and Gq. Note however that strong effects can be obtained by using peculiar symmetries 
of the rods, including the enlargment of the band gaps ||16|| . Let us now denote T (Gy) and 
r (Gq) the groups of operators associated to Gy and Go respectively fTTH - Both operators A 
and d u commute with Y (Gy) and Y (Gq). However, due to the function sq, the propagation 
equation is only invariant under the intersection group Y (Gy) nr (Gq). This simple remark 
is a crucial point in understanding the invariant properties of finite crystals and it leads to 
extremely important effects in practical applications: it is a clue to controlling the directions 
of light propagation in the structure. Indeed, due the boundary of Q, the degree of the global 
symmetry group of the device is reduced and consequently, from selection rules, the number 
of forbidden directions increases. 

In order to make this reasoning more explicit, we present two numerical experiments 
obtained through a rigorous modal theory of diffraction [0,0] that takes into account all 



of the multiple scattering between rods (in case of an infinite lattice this is exactly the KKR 
method), moreover it has been succesfully compared with experiments [[BJ. This method 



allows to deal with finite structure without periodizing tricks [20| that may lead to spurious 



phenomena As a lattice group Gy, we choose the diedral group D 6 |22|], so that the 
lattice has a hexagonal symmetry. The distance between neighbouring rods is denoted by 
d. In order to create a defect mode in that structure, we open a microcavity at the centre 



of the crystal by removing a rod (see pj,P JT^j23|j2^] for studies of the properties of defects 



in photonic crystals). A defect mode appears within the first gap at a complex wavelength 
X/d = 2.264 + 0.06z (using a harmonic time-dependence of e~' tuJt ). Such a structure can be 
used as a resonator coupled to waveguides fi~4}| . 

In the first experiment, we choose the same group for the boundary as that of the lattice 
(fig. 1) (i.e. Gq = Dq). In that case the propagation equation is completely invariant under 
r (D 6 ). We plot the map of the Poynting vector modulus associated to the defect mode (fig.l) 
and the radiation pattern is given in fig. 2. Clearly the field shows a hexagonal symmetry, 
which is obvious from the invariance of the d'Alembert equation. However, when designing 
light emitting devices, one whishes to control the direction of light emission. In this example, 
there are too many directions of emissions: such a device is to be coupled, for instance to 
waveguides, and to get a good transmission ratio, one needs to concentrate the field in a 
few useful directions. As it has been stated above, the number of authorized directions can 
be reduced by reducing the global symmetry group D§ of the device. This is what we do in 
the next numerical experiment where we have changed the boundary so that it has now a 
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rectangular symmetry (Gn = {e, s x , s y , r}, where s denotes a symmetry with respect to x 
and ?/ respectively and r is a rotation of angle 7r), the device is depicted in fig. 3. In that 
particular case, the group of the boundary is contained in the group D 6 . Then the equation 
(H) is no longer invariant under T (D 6 ) but solely under T (D e ) HT (Gn) = T (Gn) which is 
strictly contained in T (Dq). All the other transformations are now forbidden. That way, we 
expect a strong reduction of the directions of propagation of the field. 

Indeed, the map of the Poynting vector of the defect mode (fig. 3) as well as the radiation 
pattern (fig. 4) shows a strong enhancement of the vertical direction by forbidding the 
transverse directions linked to the rotations and the oblique symmetries. We have designed 
a resonator that permits to couple the radiated field in up and down directions with a better 
efficiency 

It should be noted that a group theoretic analysis gives only informations on the possi- 
ble directions of emission, the actual directions on which the field concentrates cannot be 
obtained by this mean: a rigorous computation involving a finite structure is then needed. 
Nevertheless, we have demonstrated that it was possible to strongly increase the efficiency of 
resonators by simply taking into account the symmetry of the boundary of the device. This 
remark can be used rather easily in experimental situations and could lead to a dramatic 
enhancement of the output of PBG based devices. 
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Figure captions 

Figure 1: Map of the Poynting vector modulus of the defect mode. Both the lattice 
and the boundary have the same symmetry group (Gq = D 6 ). The red line represents 
the hexagonal symmetry of the boundary of the crystal. The defect mode possesses all 
transformations of the hexagonal point group. The ratio of the rod radius to the spatial 
periode is r/d = 0.15 and the optical index is n = 2.9. 

Figure 2: Radiation pattern of the defect mode for the crystal defined in figure 1. The 
radiated power is invariant by the hexagonal point group. 

Figure 3: Map of the Poynting vector modulus of the defect mode. The global symmetry 
of the crystal is given by the subgroup G^ = {e, s x , s y ,r}. The red line represents the 
rectangular symmetry of the boundary of the crystal. The defect mode is invariant under 
Go- The ratio of the fiber radius to the spatial periode is r/d = 0.15 and the optical index 
is n = 2.9. 

Figure 4: Radiation pattern of the defect mode for the crystal defined in figure 3. The 
radiated power is invariant by the subgroup Gq. 
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